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ABSTRACT 



Aims. We study the coronal magnetic field structure inside active regions and its temporal evolution. We attempt to compare the magnetic 
configuration of an active region in a very quiet period with that for the same region during a flare. 

Methods. Probably for the first time, we use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey 
(SOLIS) to model the coronal magnetic field as a sequence of nonlinear force-free equilibria. We study the active region NOAA 10960 observed 
on 2007 June 7 with three snapshots taken during a small Cl.O flare of time cadence 10 minutes and six snapshots during a quiet period. 
Results. The total magnetic energy in the active region was approximately 3 x 10^^ J. Before the flare the free magnetic energy was about 
5 % of the potential field energy. A part of this excess energy was released during the flare, producing almost a potential configuration at the 
beginning of the quiet period. 

Conclusions. During the investigated period, the coronal magnetic energy was only a few percent higher than that of the potential field and 
consequently only a small Cl.O flare occurred. This was compared with an earlier investigated active region 10540, where the free magnetic 
energy was about 60 % higher than that of the potential field producing two M-class flares. However, the free magnetic energy accumulates 
before and is released during the flare which appears to be the case for both large and small flares. 
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1. Introduction 

Methods have been developed to extrapolate the observed pho- 
tospheric magnetic field vector into the corona. Using the fact 
that the magnetic field is dominant in solar active regions 
(ARs), we are able to neglect non-magnetic forces and to as- 
sume that the coronal magnetic field is force-free. Different 
instruments provide photospheric vector magnetograph data, 
which are used as input to the extrapolation methods. These 
data have, however had a rather low time cadence. Data of high 
time cadence are required to investigate in detail, for exam- 
ple the different evolutionary stages of solar flares. Suitable 
data include magnetic field observations of the Sun provided 
by the SOLIS Vector-SpectroMagnetograph. With a time ca- 
dence of 10 minutes, the instrument is designed to mea- 
sure multiple area scans of ARs, which enables us for the first 
time to investigate the evolution of the coronal magnetic field 
energy with a high time cadence. Many existing studies deal 
with the extrapolation based on ve ctor magnetograph data. For 



Thalmann & Wiegelmann (l2008h jerformed extrapolations o ' 
Solar Flare Telescope Vector Magnetograph data. In all of these 
studies, only one snapshot was however used or, as in the last 
of the aforementioned studies, a sequence of vector magne- 
tograms with a low time cadence of one magnetogram per day. 
Therefore, an improvement is achieved by applying our extrap- 
olation technique to the high time cadence SOLIS/VSM data 
as described in the present study. 

2. Method 

2.1. Instrumentation: The SOLIS/VSM instrument 



The Vector-SpectroMagnetograph (VSM; see iJones et al 



instance, iRegnier & Priesa (120071) dealt with the photospheric 



vector magnetic field provided by the Mees Solar Observ atory 
Imaging Vector Magnetograph. IWiegeknann et al.l ( l2005h used 
spectropolarimetric data recorded with the Tenerife Infrared 
Polarimeter of the German Vacuum Tower Telescope, and 



20021) on the Syno ptic Optical Long -term Investigations of 
the Sun (SOLIS; see iKeller et aljEool has provided magnetic 
field observations of the Sun almost continuously since August 
2003. The instrument is designed to measure the magnetic field 
vector everywhere on the solar disk. Full disk vector observa- 
tions are completed at least weekly. Multiple areas scans of 
ARs have also been available since November 2006. In addi- 
tion, longitudinal magnetic field measurements in the photo- 
sphere (at the Fe i 630.15 nm and 630.25 nm spectral lines) 
and chromosphere (at the Ca n 854.2 nm spectral line) are avail- 
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Fig. 1. Longitudinal component of the SOLIS/VSM data during 
the Cl.O flare in panels (a) - (c) and after the flare in panel (d). 

able on a daily basis. Quick-look data (JPEG images and FITS 
files) of the m agnetic field vector in a nd around automatically 
selected ARs (IGeorgoulis et al.l l2008h are available online to 
the community. The Quick-look data provides es timates of the 
magn etic field strength, inclination, and azimuth dHennev et al. 
20061 and references therein), which should, because of the 
high field strength in ARs, be comparable to fully inverted data 
only differing by a few percent. The azimuth 180°-ambiguity 
is solved using the Nonpotential Mag netic Field Calculation 



method (NPFC; see lGeorgoulisll2005h . which does not intro 



duce any error in the azimuth or any other quantities. Tools for 
full Milne-Eddington inversion are being developed to provide 
more accurate magnetic data especially in weak field regions. 
In the time period around an Cl.O flare on 2007 June 7, three 
SOLIS/VSM vector magnetograms were available to use. One 
in the rising phase of the emission, one at the time when 
the flare peaked, and one in its decaying phase (at 17:10 UT, 
17:20 UT, and 17:29 UT as shown in panel (a), (b), and (c) in 
Fig.[T] respectively). All the other magnetic field measurements 
on 2007 June 7 (between 20:30 UT and 21:42 UT) allow us to 
investigate the magnetic field structure in a period of low solar 
activity. 

2.2. Numerics: Nonlinear force-free extrapolation 

The basic equations for the computation of the nonlinear force- 
free magnetic field vector B are 



(V X B) X B = 0, 
V ■ B = 0, 



(1) 
(2) 



where Eq. ([U expresses that the Lorentz force is forced to van- 
ish (as a consequence of j || B, where j is the electric cur- 
rent density) and Eq. ^ describes the absence of magnetic 
monopoles. For reviews on how to solve these equations, we 



refer the re ader for example to ISakurail (Il989h . lAmari et al 
(Il997h . and lWiegelmamil (l2008l) . 

A special form of the force-free fields are potential mag- 
netic fields which can be computed from the longitudinal pho- 
tospheric magnetic field alone and correspond to the mini- 
mum energy state for given boundary conditions. We calcu- 
late the potential fiel d with the help of a Fast-Fourier method 
(lAlissandrakisI 1 1 98 ih . In an AR, only the energy exceeding 
that of a potential field - the so-called free magnetic energy 
- can partly be transformed into kinetic energy during dynamic 
events. Therefore, nonlinear force-free (NLFF) field models are 
required for a realistic estimation of the coronal magnetic field. 
Some existing me thods for computing NLFF fields were teste d 
and compared bvlSchriiver et all (l2006l) . iMetcalf et all ( l2008h . 
and Schriiver et al. ( 20081) . Thes e works revealed t hat th e op- 
timization meth od, proposed bv IWheatland et al. I (I2OOOI) and 
implemented by Wiegeknannl ( 20041) . was a reliable and fast 
algorithm. This approach evolves the magnetic field to repro- 
duce the boundary, force-free, and divergence-free conditions 
by minimizing a volume integral of the form 



J w(x,3;,z)(B"2|(VxB)xBpH-|V-B|2) d\. 



(3) 



where V denotes the volume of the computati onal box and 



w(x, y , z) is a weighting function (for details, see WiegelmannI 
2004 . The SOLIS data were preprocessed (for details, see 



Wiegelmann et al. 2006h so that the forced photospheric 
boundary became closer to a force-free state to provide suit- 
able boundary conditions for the minimization of Eq. Q. For 
the corresponding potential and force-free magnetic field, we 
can then estimate an upper limit to the free magnetic energy 
associated with coronal currents of the form 



1 

2/Uo 



3^ , 

pot) 



(4) 



where jiQ denotes the magnetic permeability of vacuum, and 
Bpo, and B„/// represent the total energy content of the poten- 
tial and NLFF magnetic field, respectively. To estimate the un- 
certainty in the numerical result, the code was applied to the 
original SOLIS data to which random, artificial noise had been 
added in the form of a normal distribution of amplitude ap- 
proximately equal to 1 G in the longitudinal and 50 G in the 
transversal component. The chosen noise amplitudes relate to 
the sensitivity of the VSM instrument. It measures the Stokes 
V parameter far more accurately than the parameters Q and U. 
While the longitudinal field is proportional to V, the transverse 
component is derived from Q and U, which are the principal 
source of uncertainty. 

3. Results 

3.1. Flare activity of NOAA AR 10960 

The solar activity during the week of 2007 June 4 was domi- 
nated by NOAA AR 10960. An M8.9 flare occurred on June 4 
and an Ml.O flare fired off on June 9. Furthermore, 12 C-class 
flares were detected during this week and originated in this 
group or from its vicinity. The peaks in the measured solar soft 
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Table 1. Magnetic energy of the extrapolated field. 



18:00 20:00 
Start Time ( 2007 June 07 15:00 UT ) 



Fig. 2. Solar SXR flux on 2007 June 7 in the wavelength 
range of 0.1 - 0.8 nm. Vertical lines indicate the availability 
of SOLIS/VSM data. 



X-ray (SXR) flux indicated only one C 1 .0 flare on June 7, peak- 
ing at 17:20 UT, which is from interest for the present study 
(see Fig. 121). SOUS data with a high time cadence are available 
only for June 7. 

3.2. Global magnetic energy budget 

For all magnetic field configurations, we find that the energy 
of the extrapolated NLFF field exceeds that of the potential 
field (i.e., £„/// > Epot), both being approximately 10^^ J (see 
Table [TJ. This is also the case when considering the evaluated 
relative error in the energy estimation of about 0.4 % for the 
potential and 1 % for the NLFF field (i.e. Ep„, ± 0.013 x 10^^ J 
and E„iff + 0.032 x 10^^ J, respectively). The available free 
magnetic energy is always approximately 10^"^ J with a rela- 
tive error of about 14 % (i.e. Efree ± 0.026 x 10^"^ J). These 
uncertainty ranges were checked by comparing similar results 
for 3D fields during the quiet period, and once by calculat- 
ing the energy variation in the force-free fields after adding 
artificial noise to the original magnetograms. Both Epoi and 
Eniff were highest in the phase of increasing emission from 
the Cl.O flare. During the 20 minute time period of the flare 
(17:10 UT - 17:29 UT, see Fig. |3]l the magnetic energy de- 
creased by AE„,ff = 2.01 x 10^"^ J, i.e. ^ 38 % of the available 
free magnetic energy was released. Although the flare was al- 
ready declining in intensity at 17:29 UT, it still showed a SXR 
flux of above background B-level (see Fig. |2|. The next vec- 
tor magnetogram snapshot was acquired only 3 hours later at 
20:30 UT (see panel (d) of Fig. [T]) and the free magnetic en- 
ergy had decreased further, such that x 88.16 % of the origi- 
nal amount of free energy had been released. At 20:30 UT, the 
magnetic energy was only ^ 0.6 % of the total energy and con- 
sequently the magnetic field was almost potential. 
From Fig. |2] we can see that AR10960 showed only back- 
ground B-level activity (i.e. a SXR emission < 10"^ Wm"^) 
at about 18:15 UT and almost the entire free magnetic energy 
may have been released by that time. Unfortunately, no vector 
magnetograph data was available immediately after the declin- 
ing phase of the flare and we are unable to confirm this sup- 
position. For the AR studied here, the maximum excess energy 
of a NLFF field over the potential field was about 5 % dur- 
ing the investigated period. Since this so-called free energy is 
an upper limit to the available energy to drive eruptive phenom- 
ena, consequently only a small CI .0 flare was recorded. No fur- 
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Energy of the potential and NLFF field, upper limit of the 
free energy, and excess energy of the NLFF over the potential field. 




19:»0 2a-.m 
suit Time ( 2007 June 07 17:00 UT ) 



Fig. 3. Upper panel: magnetic energy of the potential (gray) 
and NLFF (black) field. Lower panel: upper limit for the 
free magnetic energy (shown on logarithmic scale). Solid and 
dashed lines represent the recorded Cl.O flare and CME, re- 
spectively. 



ther flares occurred between 20:30 UT and 2 1 :42 UT for which 
SOLIS data is available, but five C-class flares were recorded 
about 3 hours later on 2007 June 8 between 01:00 UT and 
16:00 UT. However, a significant amount of free magnetic en- 
ergy accumulated again during the quiet period after 20:30 UT 
(see Table [1] and Fig. [Sj so that the energy content of the field 
increased and became, with some fluctuation, comparable to 
that measured before the Cl.O flare. 

From a visual inspection of the magnetic field lines within the 
extrapolation volume, we recognize some changes in the mag- 
netic field structure during the Cl.O flare (see Fig.lHi. We find 
that the field lines show their highest vertical extent when the 
Cl.O flare peaked (panel (b)). Ten minutes earlier a compara- 
ble field structure was found (panel (a)), but with a lower ver- 
tical extend. Nine minutes after the flare peaked, the field line 
configuration reached its lowest vertical extent (panel (c)), and 
more field lines of the NLFF field left the extrapolation volume. 
However, at 20:30 UT, which corresponds to the field configu- 
ration of the lowest energy content, the field clearly appears to 
have restructured, reaching on average its lowest altitude. The 
return to an almost potential structure could be assigned to a 
coronal mass ejection (CME) recorded by the SoHO/LASCO 
instrument on 2007 June 7 around 17:30 UT which could have 
bodily removed magnetic helicity of the coronal field. 
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Fig. 4. Panels (a), (b), and (c) show the magnetic field configu- 
ration during the Cl.O flare. Panel (d) shows the minimum en- 
ergy configuration. Shown are field lines of the potential (gray) 
and NLFF (black) field. For improved visibility, the z-axis is 
drawn elongated. 



4. Discussion 

We have investigated the coronal magnetic field associated 
with the NOAA AR 10960 on 2007 June 7 by analyzing 
SOLIS/VSM data. Three vector magnetograms with a time ca- 
dence of X! 10 minutes were available to investigate the mag- 
netic energy content of the coronal field during a Cl.O flare, 
and six further snapshots were acquired to analyze a very quiet 
time about three hours after the flare. Before as well as after the 
small flare, the magnetic field energy was ~ 3 x 10^^ J. 
The NLFF field had a free energy of E free « L5 X lO^^ J bcforc 
the flare. As a consequence of the flare/CME, this free mag- 
netic energy reduced by almost a factor of 10 and produced an 
almost potential configuration. Six snapshots acquired within 
a time period of about 70 minutes, during a quiet period of 3 
- 4 hours after the flare, showed again an increase in the free 



magnetic energy. Since the estimated free magnetic energy re- 
mained only about 5 % of the total energy content, no large 
eruption was produced by AR 10960. 

This is clearly different from the flaring of AR 10540 observed 
on 2004 January 18 - 21, which was analyzed in a previous 
work with the help of vector magnetograph data from the Solar 
Flare Telescope in Japan of time cadence of about 1 day. In 
this AR, the free energy was Efree ~ 66 % of the total energy, 
wh ich was sufficiently hi gh to power a M6.1 flare (for details 



see lThalmann & Wiegelm ann 2008). The activity of AR 10540 
investigated earlier was significantly higher than for the data 
analyzed in the current paper, as was the total magnetic energy. 
However, despite these differences, we also found some com- 
mon features. Magnetic energy accumulates before the flare 
and a significant part of the excess energy is released during 
the flare. The high amount of free magnetic energy available in 
AR 10540 produced M-class flares, while the relatively small 
amount of free energy in AR 10960 powered only a small C- 
class flare. In both cases, all three components of the vector 
magnetogram changed during the flare, but the energy decrease 
in the NLFF field was always higher than that of the potential 
field, i.e. the energy release was more related to the change in 
the transverse magnetic field components - which correspond 
to the field aligned electric currents in the corona - than to that 
of the longitudinal component. 
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